Late Light Curves of Type la Supernovae 



O 



P.A. Milne^ 

Naval Research Laboratory, Code 7650, Washington, DC 20375 

L.-S. The, M.D. Leising 

^ . Department of Physics and Astronomy 

(N : Clemson,S.C., 29634-0978 

<d ; ABSTRACT 

O ■ 

We extend earlier efforts to determine whether the late (t>60d) light-curves of 
^ . type la SNe are better explained by the escape of positrons from the eject a or by 

^ . the complete deposition of positron kinetic energy in a trapping magnetic field. 

We refine our selection of la SNe, using those that have extensive BVRI pho- 
tometry 35 days or more after maximum light. Assuming all SNe within a given 
Ami5(B) range form a distinct sub-class, we fit a combined light-curve for all 
^ . class members with a variety of models. We improve our previous calculations of 

^ . energy deposition rates by including the transport of the comptonized electrons. 

Q . Their non-local and time-dependent energy deposition produces a correction of as 

^ , much as <0.10™ for Chandrasekhar models and <0.18™' for sub-Chandrasekhar 

. models. 

j> ! We find that applying a filter efficiency correction, derived from measured 

^ , spectra, to B, V, R, and I light-curves after day 50 can produce a consistent 

bolometric light-curve. The V band is an accurate indicator of total emission 
in the 3500A - 9700A range, with a constant fraction (~25%) appearing in the 
V band after day 50. This suggests that the V band scales with the bolometric 
luminosity, and that the deposited energy is instantaneously recycled into optical 
emission during this epoch. Varying bolometric corrections for the other bands 
are derived. We see significant evolution of the colors of SNe la between day 50 
and day 170. We suggest that this may be due to the transition from spectra 
dominated by emission lines from the radioactive nucleus, ^^Co, to those from the 
stable daughter nucleus, ^^Fe. We show that the B, V, R, and I band light-curves 
of SNe la after t>60d can be completely explained with energy deposition from 
^^Co decay photons and positrons if substantial positron escape occurs. 



Subject headings: supernovae:general-gamma rays:observations,theory 
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1. Introduction 

Type la supernovae (SNe) are an integral part of many current astrophysical inves- 
tigations. Their luminosities briefly rival entire galaxies, making them useful as distance 
indicators to large z (Riess et al. 1998; Perlmutter et al. 1999). Understanding of their con- 
tributions to galactic chemical evolution is necessary to sort out the star formation history of 
various other stellar populations in the Galaxy (Timmes et al. 1996). The escape of gamma 
rays from their ejecta may produce a significant portion of the diffuse gamma ray background 
in the 400 - 2000 keV energy range (Watanabe et al. 1999). The escape of positrons from 
their ejecta may be large enough to explain a majority of the 511 keV annihilation radiation 
measured by the CGRO/OSSE instrument (Chan & Lingenfelter 1993; Milne, The & Leising 
1999, hereafter MTL). All of these investigations rely upon a basic understanding of the SN 
la event. 

Type la SNe are widely-believed to be the thermonuclear explosion of a mass-accreting 
carbon-oxygen white dwarf (WD). The conversion of carbon and oxygen to iron-peak and 
intermediate-mass elements releases a large amount of energy initially, and subsequently the 
decay ^^Ni — ^^Co — > ^^Fe powers the supernova through gamma- and X-ray photons and 
positrons.^ The deposition and diffusion of this energy produces the optical emission that 
defines the Type la SN. The observational support for this paradigm is typified by the ability 
of a standard carbon-oxygen WD model, W7, to reproduce the hght curves and spectra of 
the SN la 1981B (Branch et al. 1985). 

During the earliest epoch after the SN explosion, the energy released in the ^^Ni — > ^^Co 
— > ^^Fe decay is efficiently trapped by the SN ejecta. The deposition of this energy creates 
energetic electrons, which thermalize (primarily via ionization and excitation of atoms) and 
recombine. The optical photons created in the recombinations diffuse to the surface, and 
escape the SN ejecta. The emerging optical spectrum is dictated by the nature of the 
diffusion of the optical photons. The expansion of the SN ejecta lowers the column density 
to the surface, decreasing both the time for optical photons to diffuse outward and the 
efficiency of trapping of gamma-rays and energetic positrons and electrons. The SN thus 
transitions from an epoch during which the energy deposition is essentially complete and 
instantaneous and the emission depends upon the diffusion of optical photons (a "diffusion- 
dominated" epoch) to an epoch during which the diffusion time-scale is negligibly short and 



^NAS/NRC Resident Research Associate. 

^Although neutrinos are also created in these decays, the lower neutrino opacity makes neutrinos an in- 
significant contributor to the energy deposition rate during the epoch of interest in this work. See Nadyozhin 
(1994) for major features of this decay chain. 
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the emerging emission depends upon the transport of the decay products (a "deposition- 
dominated" epoch). 

Observations of SNe during these two epochs probe different characteristics of the SN 
explosion. As will be described in the next section, observations of SNe la during both 
epochs have led to considerable advances in the understanding of the SN event. This work 
concentrates upon late emission. It follows a previous work (Milne, The & Leising 1999) 
which also investigated late emission, fitting model-generated energy deposition rates to 
photometry of 10 SNe la. That work treated the 50"^ -200^^ time period in an approximate 
fashion, concentrating upon the 200'' - lOOO'^ emission and will be referenced for discussions 
relating to gamma-ray. X-ray and positron transport and the yield of escaping positrons. This 
work performs similar comparisons, but concentrates more on the 50'' -200'' period, which 
spans the transition from diffusion-dominated emission to deposition-dominated emission.^ 
In addition, this paper treats the information that is available from the R and I band 
observations of SNe la, as well as from late spectra. We compare model-generated energy 
deposition rates with B,V,R, and I band photometry for a large collection of well-observed 
SNe la, before and after creating bolometric corrections from a collection of late SN la 
spectra. Through these comparisons, we address four specific questions: 1) Is there any 
order to the late light curves of SNe la? 2) Is there observational support of the suggestion 
that positrons escape the SN ejecta? 3) Do the late light curves afford any insight into the 
correct explosion scenario (s)? 4) Are the late hght curves of SNe la suitably explained by 
the decay products of the ^^Ni — > ^^Co — > ^^Fe without the inclusion of additional sources 
of energy deposition. 



2. Physics of SN la Emission 

2.1. Emission during the first 50 days 

SN light curves peak in brightness and fade on time-scales related to the lifetimes of 
^^Ni (r=8.8'') and ^^Co (t=111'^). By 100 days after the explosion, SN la emission has faded 
to less than 4% of its maximum luminosity. Largely because of the dimming of the emission, 
SNe la are better studied during the first 50 days than during later epochs. As a result, most 
investigations address issues relating to the diffusion-dominated epoch. Once thought to be 
a homogeneous class, more and better observations have revealed that intrinsic differences 
exist among SNe la (Phillips et al. 1987, Filippenko et al. 1992, 1992b, Suntzeff 1996). 



^In this work, all SNe la will be assumed to rise to peak luminosity in 18 days. 
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These differences are evident in tfie time evolution of the spectra, the shape of the hght 
curves, and in the absolute magnitude of the emission at peak luminosity. Inhomogeneity at 
the level observed has necessitated the enlargement of the paradigm from a single, standard 
model to families of models. 

The level of understanding of this inhomogeneity has progressed on many fronts over 
the last decade. Observationally, "twins" have been found for both anomalously bright, and 
anomalously faint SNe la, solidifying the existence of a large range (and perhaps a contin- 
uum) of phenomena.'^ The very existence of inhomogeneity has altered the investigations 
of explosion scenarios. Four explosion scenarios dominate the present paradigm. The first 
consists of a carbon-oxygen white dwarf near the Chandrasekhar mass accreting hydrogen 
or helium from a binary companion until it reaches a mass at which the core carbon ig- 
nites. If the resulting burning front accelerates to become a detonation in the outer layers of 
the WD, a "delayed detonation" results. If the burning front remains sub-sonic, the result 
is a "deflagration". These Chandrasekhar mass (CM) scenarios account for inhomogeneity 
by variations in the propagation of the burning front due to density and/or compositional 
differences in the progenitor carbon-oxygen WD. The second scenario consists of a lower 
mass carbon-oxygen WD accreting a helium shell, which becomes thick enough to produce 
a helium shell detonation. This, in turn, triggers central carbon ignition. In this sub- 
Chandrasekhar (SC) mass scenario, inhomogeneity is due to the different nucleosynthesis 
that results when the progenitor mass varies from 0.65 - 1.1 Mq. The third scenario merges 
two carbon-oxygen white dwarfs, with the more massive white dwarf accreting the compan- 
ion. As with the Chandrasekhar mass scenario, central carbon ignition results. In this case, 
the accreted envelope is carbon. The range of masses of these explosions have been suggested 
to vary from 1.2 - 1.8 M0, with a roughly constant nickel yield. This scenario is also referred 
to as a "double degenerate" explosion. The fourth scenario is not a thermonuclear explosion 
at all, but an accretion- induced collapse (AlC) of a white dwarf (cither CO or ONcMg). It 
has been argued that in some cases, electron capture within an accreting white dwarf may 
lead to a collapse rather than central carbon ignition (as with Types Il/Ib/Ic SNe). These 
events eject less total mass and less nickel than the other scenarios, and have been suggested 
to explain sub-luminous SNe la. Famihes of SN models have been developed within these 
scenarios (including variations of these scenarios such as pulsed delayed detonations) . Using 
these models, many authors have demonstrated that compositional and kinematic differences 
within families of these SN models can roughly simulate spectral variations in observed SNe. 
These differences are seen and simulated near peak luminosity (Mazzali et al. 1992, 1995, 



"^Li et al. 2000 concluded that only 64% of all SNe la are normally-luminous, with ~16% sub-luminous 
and ^20% super-luminous. 
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1997, Hoflich et al. 1996, Jeffery et al. 1992, Baron et al. 1996), and during the later, 
nebular epoch (Ruiz- Lapuente et al. 1992, 1996, Liu et al. 1997a, 1997b, 1998, Bowers et 
al. 1997). 

The existence of inhomogeneity affects the use of SNe la as distance indicators at the 
level currently employed. There is an apparent absolute luminosity /peak- width relationship 
inferred from distance estimates to the host galaxies of nearby SNe la, quantified variously 
as Ami5(B) (Phillips 1993), MLCS (Riess et al. 1996), and stretch (s) (Perlmutter et al. 
1997). The ability to theoretically explain the luminosity/peak- width relationship is critical 
for relating well-observed, nearby SNe to SNe at the cosmological distances at which SNe 
la are used as distance indicators. This task is difficult because calculating the transport of 
optical photons through SN ejecta is a very complicated procedure. The ejecta is constantly 
evolving in both density and temperature. A photon crossing this non-equilibrated ejecta 
is red- (or blue-) shifted relative to local matter, wreaking havoc with the line-dominated 
opacity. Further complicating matters are the uncertainties of the cross sections of many rel- 
evant interactions. Despite these difficulties, detailed studies have been performed. Hofiich 
(1995) fit the B,V,R,1 band observations of the well-observed SN 1994D with a variety of 
Chandrasekhar mass models, demonstrating a high level of discrimination between mod- 
els. Hoflich & Khokhlov (1996) then generated B-M band light curves for Chandrasekhar 
mass and sub-Chandrasekhar models, and fit B,V,R,I band observations of 26 SNe la. One 
conclusion from that study is that sub-Chandrasekhar models appear too blue at peak to 
explain observations of sub-luminous SNe la. The same conclusion was reached by Nugent 
et al. (1997). Pinto & Eastman (1996, 2000a,b,c) investigated the influence of progenitor 
mass, nickel mass, nickel distribution, explosion energy and opacity upon the bolometric 
light curves. They then concentrated upon the ability of Chandrasekhar mass models to 
reproduce the B,V,R luminosity/peak- width relation in the range 0.85 < Ami5(B) < 1.75. 
A significant achievement of those efforts is the agreement between their calculated 18^^ -20'^ 
rise-times of the B and V peak luminosity, with the observed rise-times of SNe la (Riess et 
al. 1999a,b; Aldering, Knop & Nugent 2000).^ 

Collectively, these investigations have reached the loose consensus that the Chandrasekhar 
mass explosion scenario is the favored scenario to account for the range of optical observa- 
tions. Whether any subset of SN la events occur as sub-Chandrasekhar, merger or AIC 
explosions remains unclear. In this work, we simulate light curves for models representative 



^The models shown in Hoflich (1995) and Hoflich & Khokhlov (1996) rise to V peak in IS*^ -IS"^. Recent 
simulations by Hoflich et al. (1998) show that the V peak can be delayed by ~3 days depending upon 
the C/0 ratio, but fits of these models to observed SNe (equivalent to the previous works) have not been 
published. 
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of all four scenarios, and demonstrate general results of late light curve studies. 

2.2. Emission after 50 days 

Neither Hoflich & Khokhlov (1996) nor Pinto & Eastman carried their photometric 
simulations beyond 120*^, both groups suggesting that their simulations become inadequate 
at late epochs. As described in MTL, during the time interval of interest in this work, the 
energy deposition is dominated by interactions involving the decay products of the ^^Co — > 
^^Fe decay; gamma-ray photons and positrons. The photons possess -^30 times more energy 
per decay than the positrons, but are more penetrating. This leads to the photons initially 
dominating the energy deposition but transitioning at later times to positron dominance. The 
photon transport was performed with a Monte Carlo algorithm adopting the prescription of 
Podznyakov, Sobol, & Sunyaev (1983). A detailed description of the Monte Carlo algorithm 
is given in The et al. (1990), with its apphcation to SNe la hght curves in Burrows & 
The (1990).^ The dominant interaction is Compton scattering, which produces energetic 
electrons. These secondary electrons have a mean energy typically 300 keV at 100*^. MTL 
assigned these electrons zero lifetimes, depositing the energy in-situ. This work improves on 
that by including the transport of the secondary electrons with the same algorithm employed 
for the positron transport, but using Moller scattering rather than Bhabha scattering. The 
effect of the secondary-electron transport is small but non-negligible; in the Chandrasekhar 
mass models with radial positron escape and low ionization, the light curves are fainter by 
< 0.10"*; for low-mass sub-Chandrasekhar models, the effect reaches 0.18"*. 

Positron transport depends upon the nature of the magnetic field. Three scenarios have 
been suggested to model the magnetic field (Ruiz-Lapuentc & Spruit 1998, hereafter RLS). 
The first suggests that the field is too weak to confine positrons, and positrons follow straight 
fine trajectories, with a fraction escaping the ejecta (Colgate, Petschek & Kriese 1980). The 
second suggests a stronger field that confines positrons, but with the field lines radially- 
combed by the homologous expansion. The positrons spiral along these radial field lines with 
their pitch-angles decreasing due to the field gradient (beaming), with a fraction escaping 
the ejecta (Chan & Lingenfeltcr 1993). The third scenario suggests a strong field that is 
turbulently disordered such that positrons mirror frequently with no net transport (Axelrod 
1980). Positrons may survive non-thermally at late times, but none escape. Colgate, et al. 
(1980) argued that the first two situations arc equivalent. Simulations by MTL found this 
to be approximately tiTic for most SN models. The first two scenarios will herein be referred 



The resulting gamma-ray light curves agree with Hoflich (1995) and Hoflich et al. (1996) to within 5%. 
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to as the "radial" scenario, the third will be referred to as the "trapping" scenario. The 
differing positron escape leads to differences in model energy deposition rates, which makes 
late observations of SNe la a probe of the photon and positron transport. 

The positron and secondary electron transport was performed with a Monte Carlo al- 
gorithm, as explained in MTL. The dominant energy loss mechanisms are ionization and 
excitation of bound electrons for low levels of ionization, and plasma excitation for higher 
ionization. The model-generated energy deposition rate for the sub-Chandrasekhar model, 
HED8, is shown in Figure 1. The dashed line (D) assumes instantaneous deposition of all 
decay energy, an assumption that is invalid due to gamma-ray photon escape by 30^^. The 
dotted line (G) uses the results of the gamma-ray/sccondary electron deposition only and 
assumes no deposition of positron kinetic energy. By 100*^, the deposition of positron kinetic 
energy is an important contributor to the total energy deposition rate, but the positron life- 
times are short enough that all positron curves approximate instantaneous, in-situ deposition 
of positron energy (In). At later times, densities in the ejecta are low enough to allow escape 
in the radial scenario (dark shading, R) and non-thermal survival in the trapping scenario 
(light shading, T). Both scenarios were calculated for a range of ionizations, from 1% of all 
nuclei being singly ionized, to all nuclei being triply ionized. Positron escape removes energy 
from the ejecta, leading to fainter light curves than in-situ deposition. Lower ionization 
permits more positron escape and thus is the fainter edge of the radial curve. A trapping 
field permits large non-thermal lifetimes, but no escape. This leads to storage of energy and 
a late light curve that is brighter than in-situ deposition. For this scenario, low ionization 
leads to a larger energy storage and the brighter edge of the trapping curve. ^ 



2.3. Model/Observation Compcirisons 

The ideal model/observation comparison would entail an explicit NLTE determination 
of photon and positron energy deposition rates and from these rates the generation of op- 
tical/IR spectra. These spectra would then be compared with a sequence of observed SN 
la spectra, comparing the evolution of the total flux from early to very late epochs. This 
ideal situation is far from being realized, both computationally and observationally. Falling 
short of that solution, an alternative approach has been employed in numerous works; fit- 
ting energy deposition rates to multi-band photometry or estimates of the uvoir bolometric 



^The normally- luminous sub-Chandrasekhar model shown allows a relatively large escape fraction for the 
radial scenario compared to normally-luminous Chandrasekhar mass models. Nonetheless, the basic features 
shown in Figure 1 are characteristic of all the SN la models simulated. 
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luminosity. Two issues must be addressed for these comparisons to be meaningful, photon 
diffusion delay and color evolution. Photon diffusion delay refers to the portion of the time 
delay between the energy deposition (from the scattering of gamma- and X-ray photons and 
slowing of positrons) and the production of the optical light that can be seen by an external 
observer that is due to the diffusion of optical photons out of the ejecta. This diffusion delay 
determines the shape of the SN light curve peak at early times, but decreases to negligible 
values at late times. Thus, at late epochs, energy deposition is instantaneously recycled into 
optical emission. In previous studies of late light curves of SNe la, two approaches have been 
used to account for photon diffusion delay. Studies that attempt to fit observations continu- 
ously from early to late epochs derive this time delay as a function of explosion epoch. The 
alternative approach is to study only the late epochs and assume instantaneous recycling. 
Observationally, the late epoch begins when the emitted spectrum transitions from contin- 
uum emission to (forbidden) nebular lines. Color evolution refers to the changing fraction 
of the total optical flux that a given photometric band samples due to the evolution of the 
spectrum. Again, two different approaches have been employed. The first approach cre- 
ates bolometric light curves from the photometry (applying weighting factors derived from 
spectra) and compares these light curves to model-generated energy deposition rates. ^ The 
second approach determines the epoch after which a given photometric band scales with the 
total optical fiux (the epoch is determined from spectra), and compares this band photometry 
to model-generated energy deposition rates. 

Two groups that carried their investigations to late epochs have treated photon diffusion 
in their comparisons. Colgate, Petschek and Kriese (1980) fit Monte Carlo model simulations 
of the energy deposition rate to B band data from SNe 1937C and 1972E. That work assumed 
a grey opacity for both photon and positron transport and a single-zone model. They 
concluded that positron escape was required to explain the late light curves. Cappcllaro 
ct al. (1998) (hereafter CAPP) and Salvo ct al. (2001) transported gamma-ray photons, 
positrons and optical photons through W7-Uke models (scaled to various masses) with a 
Monte Carlo code and fit the energy deposition rates to V band data from SNe 1991T, 
1994D, 1992A, 1993L, 1996X and iggibg.'^ That work assumed a grey opacity for photons 
and varied the positron opacity to fit the late light curves. They concluded that positron 
escape was suggested for some, but not all SNe. 

RLS transported gamma-ray photons and positrons through Chandrasekhar mass and 



^Contardo et al. (2000) discuss the issues related to generating bolometric light curves from band pho- 
tometry. 

^Both Colgate et al. (1980) and CAPP transported the positrons with routines that were produced for 
photon transport. The positrons were given a speed c. 



-9- 



sub-Chandrasekhar models without treating photon diffusion delays, and fit the energy de- 
position rates to bolometric light curves for the SNe 1972E, 1992A and 1991bg.^° That 
work concentrated upon very late times (> 100*^) at which time the photon diffusion delay 
time-scale is assumed to be short. They concluded that normally-luminous SNe la transi- 
tion from a positron-trapping phase to a positron-escape phase after 450^^, and can be fit 
with Chandrasekhar mass models. They further concluded that sub-Chandrasekhar mod- 
els best describe the sub-luminous SNe la and feature positron escape as early as ISO''. 
MTL transported gamma-ray photons and positrons through Chandrasekhar mass and sub- 
Chandrasekhar models, and fit the energy deposition rates to bolometric, V band and B 
band light curves for 10 SNe, including the SNe studied by the previously mentioned groups. 
That work also concentrated on very late times and did not treat photon diffusion delays. 
They concluded that both Chandrasekhar mass and sub-Chandrasekhar models can explain 
normally- and super- luminous SNe la, with positron escape suggested to be consistent with 
seven of the eight SNe. They further concluded that none of the models tested could explain 
the sub-luminous SN 199 Ibg, a result in agreement with CAPP, but in disagreement with 
RLS. 

This work docs not simulate photon diffusion delays, comparing model-generated en- 
ergy deposition rates to multi-band photometry as though the delays were negligible. This 
approach uses the comparisons to determine the onset of instantaneous recycling. Similarly, 
the initial comparisons of the model- generated energy deposition rates to each of the B,V,R 
and I band data-sets assume no color evolution. The differences between these comparisons 
and the later comparisons that incorporate spectrally-derived bolometric corrections afford 
a measure of the importance of color evolution at late epochs. 



3. SN la Observations 

To perform model-generated energy deposition rate/SN light curve comparisons, a large 
collection of SNe la observations have been compiled. The principle sources of B,V,R,I band 
photometry are the 29 SNe la observed by Hamuy et al. (1996), and the 22 SNe la observed 
by Riess et al. (1999). MTL fit the 10 SNe la best observed to very late times. Of those 
ten SNc, only six are included in this study (as explained below). These SNe have been 
supplemented with observations of SN 1995D by Sadakane et al. (1996), SN 1996X by Salvo 
et al. (2001), SN 1997cn by Turatto et al. (1998), SN 1998bu by Jha et al. (1999), Suntzeff 



^°The bolometric light curves were generated from photometry for SNe 1992A and 1991bg. For SN 1972E, 
the bolometric light curve was generated from a series of optical spectra. 
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et al. (1999), Garnavich et al. (2000), Cappellaro et al. (2001), and SN 1998de by Modjaz 
et al. (2000). The SNe la have been divided into luminosity sub-classes (normally-, sub- 
and super- luminous). Hereafter these categories will be denoted in figures by N, sb 8z SP. 
The categorization is based upon the Ami5(B) (Phillips 1993) values of the SNe. Numerous 
empirical relations have been derived for SNe la, plotting a given variable versus the Ami5(B) 
parameter. Shown in Figure 2 are the maximum B,V and 1 band absolute magnitudes for 21 
of the 22 SNe la analyzed in this study combined with 18 others from Hamuy et al. 1996.^^ 
Early studies assumed a hnear relationship between these variables for all SNe la. A larger 
samphng has suggested that all these relations deviate from linearity at Ami5(B) ~ 1.6. 
It is unclear whether this deviation is best explained by a continuous quadratic function 
applicable to all SNe la (Phillips 1999), or by two different linear relations, one for each sub- 
class of SNe la. Regardless of the interpretation of this break, the normally- /sub- luminous 
cut has been set at 1.60. By that definition, SN 1993H is considered sub-luminous and SN 
1992A is considered normally-luminous. The normally-/supcr-luminous cut (Ami5(B)=0.95) 
was chosen such that the well-observed SN 1991T (which showed spectral deviations from 
normal SNe la) was considered super-luminous. It is important to note that the absence 
of a template SN with a wider peak than SN 1991T hinders the accurate estimation of 
Ami5(B) for super- luminous SNe. As will be shown, the exact location of the normally- 
/super-luminous cut is unimportant for this work, but the normally-/sub-luminous cut does 
differentiate two distinct late light curve sub-classes. 

Prom 64 observed SNe la, 22 were included in this study. To be included, a SN 
had to be observed at least once after SS*^ with at least one observation within ±20 of 
the normalization epoch (explained below), chosen to be 65*^ post-explosion. Addition- 
ally, the SN had to be first observed no later than one week post-B maximum (a re- 
quirement relaxed for the SN 1992K due to the undersampling of the sb sub-class). The 
SNe la which meet this criteria arc; N:90N,90O,92A,92al,93ag,94D,95D,95E,95aL96X,98bu, 
SP:91T,91ag,92bc,94ae,95bd, sb:86G,91bg,92K,93H,97cn,98de. Four SNe la from MTL were 
not included in this study. SNe 1937C and 1972E were not well-observed with multi-band 
photometry and were excluded. SN 1993L was discovered more than a week post-peak and 
was excluded. The late photometry of SN 1989B is being re-analyzed and thus that SN was 
excluded. All three sub-classes are better sampled from SO*^ -150'^ than after ISO''. 

The phenomenon studied in this work is the evolution of the light curves at late epochs. 



^^This figure only included the SNe with well-defined peak shapes. 

^^The SN has been identified as a light echo candidate (Bofii et al. 1999). As published, the late light 
curve of SN 1989B is more than 1™ brighter than the suggested models and is the singular exception to the 
trends seen in the other SNe. 
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For this reason, a relative magnitude scale has been employed rather than an absolute or 
apparent magnitude scale. All SNe have been normalized to have zero magnitude at 65*^ 
post-explosion (assuming an IS*^ rise-time for all sub-classes). The initial motivation for the 
65'^ normalization date was based upon the suggestion in Pinto & Eastman (2000) that the 
photon diffusion time-scale becomes negligibly short at roughly that time. As the next section 
shows, a high level of homogeneity exists within the multi-band photometry of each sub-class 
after 65*^. Additionally, the greater number of 45'* -85'* observations relative to later epochs 
makes the normalization algorithm perform better at early times. These two factors led us 
to retain the 65*^ normalization date for all data-sets, although the model-generated energy 
deposition rates do not fit the band photometry at GS'^ for any band/sub-class combination. 
The normalization has been performed by linear interpolation of all data between 45*^ and 
85^^. For two SNe, only a single observation existed in that interval. For those SNe, the 
latest pre-45'^ observation was used to perform the interpolation. 

By treating all SNe within a given Ami5(B) range as a single object, we are assuming 
homogeneity within each luminosity sub-class. The amount of scatter about a template gives 
some measure of the quality of that assumption, but numerous sources of systematic errors 
exist to cloud that interpretation. Each SN light curve is measured relative to the host 
galaxy's background light, which is unique to that SN. If the SN is near the detectability 
limit for that observation, the infiuence of the background subtraction increases in impor- 
tance. The irregularity of the sampling means that highly uncertain observations are mixed 
throughout every epoch of the late hght curves. Similarly, different telescopes, filters and 
detectors were used to produce this data-set, introducing additional systematic differences 
between individual data-points. As was shown by Suntzeff (1999), observations of the same 
SN by the same observer at the same site with the same class telescope can yield different 
photometry. The differences were found to be in excess of the quoted uncertainties. Suntzeff 
showed that in this optimal case, the "instrumental uncertainty" was on the order of 0.06™ 
for the V band during the GO'' -75'^ epoch. For the case of this mixed data-set, the uncertain- 
ties are likely to be larger. Errors in this work were also introduced by the fitting algorithm, 
which fit the individual light curves to the estimated GS'^ magnitude. 

All errors discussed above have been related to estimates of the relative magnitude. As 
SNe are always detected at some time after the explosion, the explosion date is also somewhat 
uncertain. Two factors may infiuence uncertainties in the explosion date; 1) uncertainty of 
a sub-class rise-time to peak, and 2) uncertainty of the peak date. For this work, all SNe 
were assumed to rise to peak B magnitude in 18'* regardless of the luminosity sub-class. To 
minimize the effect of peak date uncertainty, only SNe first observed no later than a week 
after maximum light were used (with the exception of SN 1992K). The determination of the 
peak date depends upon template fitting, so peak date uncertainties are only relative to the 
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template used. SNe discovered before tB(max) have peak date uncertainties typically of il*^. 
SNe discovered post-peak have larger uncertainties. Nonetheless, as the time-scales involved 
in this work are on the order of tens of days, uncertainties on the order of a few days will 
not significantly alter these results. 

4. Results 
4.1. V Band 

The 22 SNe included in this study are shown fit (by eye) with the dclayed-detonation, 
Chandrasekhar mass model DD23C in Figure 3. All 22 SNe have been normalized at 65'^. 
Three tendencies are apparent in the data. First, within each sub-class the data show re- 
markable homogeneity. The sub-classes are DEFINED by differences in the early light curve 
shape, but after 65*^, there is no evidence of variations. Second, the normally-luminous and 
super-luminous data appear to evolve similarly, while the sub-luminous sub-class continues 
to fall more steeply (recall that a steep decline from peak DEFINES this sub-class). Third, 
the separation between the normally-/super- luminous and the sub-luminous sub-classes far 
exceeds the scatter within the sub-classes. This argues against a continuous transition be- 
tween the sub-luminous and normally-luminous sub-classes. Unfortunately, the sub-luminous 
sub-class is under-sampled, so many more sub-luminous SNe need to be observed before ad- 
dressing whether there exists an absolute separation of sub-classes. 

After ^80'' the normally-luminous and super-luminous sub-classes appear to be explain- 
able by this model, assuming radial escape. By contrast, the sub-himinous sub-class cannot 
be explained by either magnetic field scenario. To improve the visualization of the 40^^ -120*^ 
normally-luminous and super-luminous data, in the inset of Figure 3 we show the data as 
residuals relative to the model curve. The similarities between the normally-luminous and 
super- luminous data-sets is apparent. Also apparent in the inset are the failure of the model 
light curves to fit the V band data in detail until after ^120^^. These results are not specific 
to the delayed-detonation, Chandrasekhar mass model, DD23C. Shown in Figure 4 are six 
models suggested to explain normally-luminous and super-luminous SNe, shown in the "delta 
magnitude" format. The models are defined in Table 1. All model light curves show the 
same fundamental structure, with the models fitting the data after '^170'^ and suggesting 
positron escape. The Chandrasekhar mass models fit the normally-luminous data better 



^^The delta magnitude format shows the residuals of the data and model-generated light curves to the 
instantaneous deposition approximation, in units of magnitude. 
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than does the sub-Chandrasekhar mass model, but the improvement is modest. As was 
discussed in MTL, the super-luminous model light curves are similar, both to other super- 
luminous SN models and to normally-luminous Chandrasekhar mass model light curves. We 
assert that the 16 normally- and super-luminous SNc do not need to be differentiated into 
distinct sub-classes after 60^^. We further assert that if the V band emission scales with the 
bolometric emission to the level of a few tenths of a magnitude, then the late emission from 
all 16 SNe can be suitably explained with a single light curve featuring positron escape, and 
that Chandrasekhar mass, sub-Chandrasekhar mass, and merger explosions can all explain 
these data. 

It is apparent from Figure 3 that the sub-luminous sub-class differs considerably from the 
normally-/super-luminous sub-classes. MTL concluded that none of the models they tested 
could explain the 60*^ -560'^ evolution of the light curves. One interpretation of that result 
is a rejection of both the Chandrasekhar mass and the sub-Chandrasekhar mass models 
tested. A different interpretation would be the failure of the assumption that the energy 
deposition rate scales with the V band during that epoch. Figure 5 shows five sub-luminous 
models fit to the sub-luminous sub-class data. The models are normalized to fit the data at 
a later epoch (~170''). It is clear from the figure that none of the models can explain the 
V band data before 150''.^'' After 150'^, the V band data can be suitably explained by all 
models if positron escape is assumed. This latter interpretation implies that the energy 
deposition rate scaling with the V band occurs later for the sub-luminous sub-class than for 
the other sub-classes. For this interpretation to be correct, there must be an explanation 
which is consistent with the spectral observations, which suggest that sub-luminous SNe la 
enter the nebular phase earlier than normally-luminous or super-luminous SNe la (Mazzali 
et al. 1997). 

4.2. BVRI Bands 

In the previous section, we fit model-generated energy deposition rates to V band data. 
That procedure relies upon the assertion that the V band traces the bolometric luminosity. 



^''The fifth panel (lowest) is an accretion-induced collapse model, which ejects 0.2 Mq of material. The 
model was created by scaling ONeMg to 0.2 Mq. This model is included upon the suggestion of Fryer et 
al. (1996,1999) that AlC's may explain sub-luminous SNe la. Even at this extremely low mass, the early 
escape of positrons is not enough to reproduce the steepness of the 65''' -170'' light curve. 

^^It is important to note that the range of ionizations allowed lead to fairly thick curves for low-mass 
models, which exaggerates the visual impression of their ability to explain the data. A treatment of the level 
of ionization would reduce the range to a single curve. 
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For the normally-/super-luminous data-set, the data are approximately fit with the models, 
suggesting a validation of that assertion. By contrast, for the sub-luminous data-set the 
models could not fit the data until a later epoch. All SNe la used in this study were also 
observed in the B, R and 1 bands. By independently fitting these data-sets to the model- 
generated energy deposition rates, we further explore the nature of the late-time emission 
from SNe la. 

Shown in Figure 6 are the BVRl data-sets for SNe la fit to the model, W7. Within each 
sub-class the data has been normalized to 65^^. Each sub-class has then been independently 
fit to the model light curves. For all four photometric bands, the normally-luminous and 
super-luminous light curves are similar. The B band data for the normally-luminous and 
super- luminous SNe la deviate considerably from the model before 100*^, but by ~170'^ they 
roughly scale with the energy deposition rate. The B band data for the sub-luminous sub- 
class differ from the normally-luminous and super-luminous sub-classes early, but appear 
similar at late times. Perhaps the most striking results are in the R and I bands, where the 
three sub-classes are quite similar. We argue that the differences between the late evolution 
of the B and V band data do not carry over to the longer wavelengths. The model fits initially 
deviate from the R and 1 band data for all three sub-classes, but then converge upon the 
model curve by ~170'^. The convergence in the I band is less certain due to scatter. 

It is interesting to note that the onset of the energy deposition scaling with the pho- 
tometry {^170'^) occurs later than the transition from photospheric continuum emission to 
nebular emission (^60^^) for all three SN la sub-classes and all four photometric bands. To 
determine whether this delay is explainable as a consequence of color evolution, we produce 
spectral templates for normally-/super-luminous and sub- luminous SNe la. The templates 
have been created from a collection of SN la spectra; 1994D (Patat et al. 1996; Turatto 
2000; Fihppenko 1997), 1991T (Fihppenko et al. 1992), 1996X (Salvo et al. 2001), 1991bg 
(Turatto et al. 1996; Ruiz -Lapuente et al. 1993), 1987L (Ruiz -Lapuente et al. 1993), 
1994ae (Bowers et al. 1997), 1998de (Modjaz et al. 2000), 1998bu (Jha et al. 1999; Math- 
eson 2000), 1995D (Matheson 2000), 1981B (Branch et al. 1983), 1984A (Branch 1987), 
1991F (Gomez & Lopez 1995), 1992 A (Kirshner et al. 1992). We include all emission within 
the 3500A-9700A wavelength range, and ignore all emission outside of this range. Spectra 
that do not span that range have been linearly interpolated from earlier and later spectra. 
The spectral sequences for a sampling of each sub-class are shown in Figures 7 and 8. From 
these templates, we have determined the evolution of the "filter efficiency" , or the fraction 
of the energy emission in the 3500A-9700A wavelength range that is detected by each band. 



^^The interpolation was performed in two steps. The shape of the interpolated portion was determined 
from the two complete spectra adjacent in time. That shape was then spliced to the existing portion. 
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The filter efficiencies for normally-luminous and super-luminous SNe la, shown in Figure 9, 
match the shape of the residuals of the data-model fits shown in Figure 6 for the normally- 
/super-luminous data-set. The V band filter efficiency remains nearly constant, in agreement 
with previous assertions that the V band scales with the bolometric luminosity. ^"^ Shown 
in Figure 10 are the B,V,R and 1 band data fit to the model W7, after deriving bolometric 
corrections from these filter efficiencies. The data for all four bands are suitably explained by 
the energy deposition rates of the model if positron escape is assumed. Before ~60'', the B,V 
and I data vary dramatically from the energy deposition rates. This would be consistent with 
the decrease to negligible values of the photon diffusion time-scale, as suggested by Pinto 
& Eastman (2000a). The correction lessened, but did not completely remove the residuals 
from the B and V band data during the 60*^ -170*^ epoch. It is unclear whether the residuals 
remain due to a modest failure of this relatively crude bolometric correction calculation, or 
rather they remain because this treatment has ignored some physical phenomenon. 

The application of these corrections essentially creates bolometric light curves from each 
photometric band. A single bolometric light curve for normally- and super-luminous SNe la 
was generated by combining the four curves shown in Figure 10 after weighting each curve 
according to the photometry errors. The resulting bolometric light curve can be compared 
with model-generated light curves. These comparisons, shown in Figure 11, exhibit the same 
tendencies that are seen in the V band; the energy deposition rates for all four models are 
consistent with the data, and positron escape is suggested. These curves differ slightly from 
the bolometric fight curve derived by Contardo et al. (2000) from much of the same data. 
That study was primarily concerned with emission near the luminosity peak and assumed 
constant bolometric corrections at late times (> 130*^). As we have shown that some evolution 
occurs during the 60^ - 200^^ time period, we assert that this bolometric light curve is more 
accurate at late times than the light curve from Contardo ct al. (2000). Although bolometric 
light curves include information from four photometric bands rather than simply from the V 
band, the V band is the best observed. Thus, we present the bolometric light curve primarily 
to demonstrate that conclusions derived from studies of the V band shape are valid. 

There are few nebular spectra for sub-luminous SNe la, so for this sub-class, the filter 
efficiencies are very crude and are undersampled after 150*^ (Figure 12). This is problematic, 
because as shown in Figure 6, the 100*^ -200^ time-span is a critical transitional epoch. 
Nonetheless, as shown in Figure 13, the B,V,R and I band data are also fairly sparse and 
can be roughly explained with the model PDD54, after correcting for filter efficiency. All four 
bands fit the energy deposition rates better with the correction, as seen through comparison 



Contardo et al. (2000) arrived at a similar result. 
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with Figure 6. There are large residuals in the V and I bands before ^SO''. This is hkely due 
to photon diffusion delays. It is evident in Figure 13 that the light curves of sub-luminous 
SNe la are consistent with model-generated energy deposition rates which allow positron 
escape. More spectral and photometric data of sub-luminous SNe la need to be collected to 
further the understanding of this epoch. Due to the crudeness of the bolometric corrections, 
a bolometric light curve was not generated for the sub-luminous sub-class. 

The success of the bolometric corrections in reducing the residuals during the 60°' -170°' 
time-span suggests the existence of an epoch during which optical emission is instantaneously 
recycled, but the SN color continues to evolve. A possible explanation for this color evolution 
is the decay of ^^Co to ^^Fe. It is during the 65"^ -170*^ epoch that the daughter becomes the 
dominant species in that decay. As shown by both Liu et al. (1997a,c) and Bowers et al. 
(1998), cobalt emission dominates the wavelength ranges observed with the R and I bands. 
As ^^Co decays to ^^Fe, these bands would fade faster than the energy deposition rate due to 
color evolution until late-times at which time most cobalt has decayed and iron and stable 
nickel dominate the spectrum in that wavelength range. In the V band wavelength range, 
both Co and Fe lines are present and the bolometric correction would change only slightly 
during the transition. The B band may brighten due to the increased emission from the 
many [Fell] lines from 446nm -456nm. This explanation can be tested through comparisons 
of model-generated spectra with observations made during this epoch. 



5. Discussion 

Two tasks comprised this work. The first task was the compilation of a data-set of 
B,V,R,I photometry of type la SNe divided into sub-classes, and development of bolometric 
correction factors for each band. The second task was the generation of late-time light curves 
from SN la models and the fitting of these models to the compiled data-set. Certain findings 
can be derived from this study independent of the model-fitting, and are worthy of mention. 
We have demonstrated that there is order to the late light curves. The late B and V band 
light curves of sixteen normally- and super-luminous SNe follow a similar evolution, while 
those of six sub- luminous SNe follow a different evolution. By contrast, the late R and I 
band light curves of all 22 SNe la follow a similar evolution. The normally-/super- luminous 
sub-classes appear virtually identical although they span more than half of the Ami5(B) 
range of SNe in this study. The B and V band light curves of the sub-luminous SNe appear 
to require a distinct sub-class. Whether there are in fact SNe la that connect them to 
normally-luminous objects is an important question. This question can only be answered by 
observational campaigns that focus upon monitoring SNe that span the range of Ami5(B) 
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values. 

Comparisons of model-generated energy deposition rates with photometry in the B,V,R, 
and I bands reveal three epochs. There is an early epoch (< 60''), where energy depo- 
sition/photometry comparisons are invalid presumably due to the time delay between the 
deposition of energy and the emergence of the resultant optical emission. That epoch is 
followed by an intermediate epoch (60'^ -170'^) where the bolomctric light curve is fit by the 
energy deposition rate, but color evolution must be addressed for individual photometric 
bands to be fit by the energy deposition rate. This epoch is followed by a late epoch (> 
170°') where the bolomctric corrections arc roughly constant and each individual photomet- 
ric band can be fit by the energy deposition rate if positron escape is allowed for. These 
epochs exist for both the normally-/super-luminous sub-class, and the sub-luminous sub- 
class. The bolomctric corrections during the intermediate epoch are better defined for the 
normally-/super-luminous SNe la than for the sub-luminous SNe la, the latter corrections 
are preliminary. 

Collectively, the hght curves of SNe la after 60*^ suggest that the interactions of the 
products of the ^^Co — > ^^Fe decay with the ejecta can explain the energy deposition without 
any additional energy deposition source. Positrons are seen to escape the ejecta in quantity 
for all SNe la during the late epoch. These findings hold for Chandrasekhar mass, sub- 
Chandrasekhar mass, merger and AlC models equally. The model independence of these 
results docs not imply that the late emission from SNe la cannot probe the progenitor. The 
challenge for NLTE radiation transport calculations is to reproduce the spectra during this 
epoch. 

There are two ramifications of positron escape from SN la ejecta that warrant further 
discussion. First, nebular spectra have been used to estimate the ^^Ni production in SNe 
la (Ruiz-Lapuentc & Filippcnko 1996, Bowers et al. 1997) and to differentiate between 
Chandrasekhar mass and sub- Chandrasekhar mass models (Liu et al. 1997a, 1997b, 1998). In 
all of these studies, instantaneous positron energy loss was assumed. By 300*^, the energy 
deposition is dominated by positron slowing, and a substantial fraction of those positrons 
are expected to escape the ejecta in the radial scenario. Even positrons that do not escape 
the ejecta diffuse from the location of their creation. In light of this and other recent studies, 
nebular spectrum studies should be calculated with realistic positron transport. Second, 
positrons that escape the SN la ejecta are thought to survive non-thermally on time-scales of 
10^ years, or more (Chan & Lingenfelter 1993, Guessom, Ramaty & Lingcnfclter 1991). The 
collective positron contributions from SNe la, as inferred from extra-galactic SN rates, are 
sufficient to generate a large fraction of the 511 keV positron annihilation radiation observed 
by the CGRO/OSSE, SMM and TGRS gamma-ray telescopes (Milne et al. 1999,2000). 
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The degree of dominance is enough that accurate characterization of positron annihilation 
radiation will trace recent SN la activity in the Galaxy. 

Light echoes have been detected from two SNe la, 1991T and 1998bu. In addition, BofR 
et al. (1998) have identified SN 1989B as a hght echo candidate. Of all SNe la observed 
after 300^, these three are the only SNe to remain significantly brighter than those in this 
work. As light echoes are currently being used as geometric distance indicators, this work 
suggests that the templates shown in this work may be used to locate light echoes (Sparks 
et al. 1999). 

There remain incomplete aspects to this study. More SNe la need to be monitored to re- 
duce gaps in observations and determine the late behavior of more SNe la along the Ami5(B) 
sequence. The exact nature of the transition from the onset of instantaneous recycling to 
the cessation of color evolution must be better observed spectrally and in each photometric 
band, particularly with sub-luminous SNe la. These unexplained aspects underscore the fact 
that much is to be learned about the physics of Type la SNe from late observations. 

We thank P. Hofiich for ongoing access to SN models and for discussion relating to 
radiation transport through SN ejecta. We thank M. Turratto and S. Jha for access to 
SN photometry and spectra. We also thank T. Matheson, R. Lopez, and P. Meikle for SN 
spectra. 
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Fig. 1. — A model generated bolometric light curve for the model HED8. The dashed 
hne (D) assumes instantaneous deposition of all decay energy. The dotted line (G) uses 
the results of the gamma-ray energy deposition only and assumes no deposition of positron 
energy. Between these two boundaries are the results of the gamma-ray energy deposition 
coupled with instantaneous positron deposition (thick line, In) and the range of curves for 
a radial field geometry (dark shading, R) and for a trapping geometry (light shading, T) as 
the electron ionization fraction varies from 0.01 < Xe ^ 3- 
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Fig. 2. — The absolute B,V, & I magnitudes of SNe la as a function of the Ami5(B) param- 
eter. All values taken from Hamuy et al. 1996, except SNe 1997cn (Turatto et al. 1998) and 
1998de (Modjaz et al. 2000). The normally-luminous, super-himinous and sub- luminous 
SNe la are shown as filled circles, open circles, and crosses, respectively, while other SNe not 
used in this study are shown as filled squares. All bands' absolute magnitudes show a break 
around Ami5(B)=1.6. 
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Fig. 3. — The V band light curves of 22 SNe la fit with the energy deposition rates 
of the CM model, DD23C. Normally-luminous SNe la (filled circles) and super-luminous 
SNe la (open circles) appear to follow similar evolutions at late times, while sub-luminous 
SNe la (crosses) follow a different evolution. The data is better fit at late times by 
radial escape of positrons (R: light shading) than by positron trapping (T: dark shad- 
ing). Inset: detailed 30*^ -120^ evolution of the N &; SP sub-classes (see text). SNe 
1990O,1990T,1991ag,1992K,1992al,1992bc,1993H,1993ag from Hamuy et al. 1996. SNe 
1994ae,1995D,1995E,1995al,1995bd from Riess et al. 1997. Other SNe used, SN 1992A 
(Suntzeff 1996), 1990N (Lira 1998), 1991T (Lira 1998; Schmidt et al. 1994), 1991bg (Lei- 
bundgut et al. 1993; Turratto et al. 1996; Filippenko et al. 1992), SN 1994D (Patat et al. 
1996, Tanvir 1997, CAPP, Cappellaro 1998), SN 1986G (Phillips et al. 1987, 1998, Cristiani 
1992), SN 1997cn (Turatto et al. 1998), & SN 1998de (Modjaz et al. 2000), SN 1996X (Salvo 
et al. 2001). SN 1995D also from Sadakane et al. 1996. 
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Fig. 4. — The V band light curves of 16 normally-luminous SNe la (upper four panels-filled 
circles) and super-luminous SNe la (lower two panels-open circles) fit with model-generated 
energy deposition rates from six SN la models, shown in the delta magnitude format. The 
four N and two SP models have been suggested to explain normally- and super-luminous 
SNe, respectively. The data is approximated by all six models after ITC* if radial escape of 
positrons is assumed (lower band shaded curves). The models' light-curves with positrons 
trapped (upper band shaded curves) cannot give a good fit to the late-time data. 
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Fig. 5. — The V band light curves of 6 sub-luminous SNe la fit with model-generated energy 
deposition rates from five SN la models. Models and data have been transformed to the delta 
magnitude format. The data is consistent with all five models after ITO'' if radial escape of 
positrons is assumed. 
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Fig. 6. — The BVRI band light curves of 22 SNe la fit with model-generated energy de- 
position rates from the CM model, W7. Models and data have been transformed to the 
delta magnitude format. The symbols of the data points are as in Fig. 4. The sb luminosity 
sub-class follows a different early evolution than the N & SP subclasses in the B & V bands, 
but all three sub-classes are similar in the R & I bands. For all SNe la, all four bands are 
approximated by the model-fit after 170'* if positron escape is assumed. 
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Fig. 7. — A sequence of spectra of normally- and super-luminous SNe. All spectra have 
been offset by an arbitrary constant. Incomplete spectra are denoted with an asterisk (*). 
It is clear from these spectra that the emission within this wavelength range experiences a 
blueward shift. The transmission efficiencies for the B,V,R, and I filters are shown above 
the spectra for reference. 
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Fig. 8. — A sequence of spectra of sub-luminous SNe. As in Figure 7, all spectra have been 
offset by an arbitrary constant, and incomplete spectra are denoted with an asterisk (*). 
This spectral sequence differs considerably from Figure 7, most notably due to the peak at 
7300A, and the persistence of emission red-ward of 8OOOA. The transmission efficiencies for 
the B,V,R, and I filters are shown above the spectra for reference. 
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Fig. 9. — The evolution of the fraction of the emission (3500A-9700A) observable by the 
B,V,R and I bands for N and SP SNe la. Filled diamonds show spectra that spanned 
the wavelength range, open diamonds are incomplete spectra for which interpolation was 
required. 
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Fig. 10. — Color-corrected BVRI band light curves of 16 N and SP SNe la fit with model- 
generated energy deposition rates from the CM model, W7. As in Figure 6, the models and 
data have been transformed to the delta magnitude format. 
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Fig. 11. — Model-generated energy deposition rates for four SN la models compared with 
calculated bolometric light curves for 16 normally- and super-luminous SNe la. 
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Fig. 12. — The evolution of the fraction of the emission (3500A-9700A) observable by the 
B,V,R and I bands for sb SNe la. As in Figure 8, Filled diamonds show spectra that spanned 
the wavelength range, open diamonds are incomplete spectra for which interpolation was 
required. 
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Fig. 13. — Color-corrected BVRI band light curves of 6 sb SNe la fit with model-generated 
energy deposition rates from the CM model, PDD54. As in Figures 6 & 9, the models and 
data have been transformed to the delta magnitude format. 



